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1. Introduction 

Supemovae (SNe) exploding in a dense circumstellar medium (CSM) are hypothesized to 
effectively accelerate cosmic rays (CR) in interactions of the SNe ejecta with the CSM. [1, 2] 
and [3] showed that a collisionless shock is necessarily formed after the shock breakout if the 
SN progenitor is surrounded by an optically thick CSM. Cosmic rays might be accelerated to 
PeV energies and contribute to the knee structure in the cosmic-ray spectrum [2, 4]. pp and py 
interactions produce pions, which in the neutral case decay to gamma rays and in the charged case 
produce neutrinos in the decay chain. Interaction of the ejecta with the CSM can convert a large 
fraction of the kinetic energy to radiation causing this type of SNe to be more luminous than other 
core-collapse SNe [5, 6]. Figure 1 (taken from [2]) shows a schematic picture of a SN exploding 
in dense circumstellar medium. 

Several candidates for such interaction-powered SNe were found [7, 8, 9, 10] and some su- 
perluminous supemovae were suggested to be powered by interactions [11, 12]. Observationally, 
Type Iln SNe (“n” for narrow) are the best candidates to be interacting with a dense CSM. They 
are spectrally characterized by the presence of strong narrow emission lines, mostly Ha, originat¬ 
ing from surrounding H II regions. A slow spectral evolution indicates the presence of a dense 
CSM [13, 14]. They are often accompanied by precursor mass-ejection events, which could be 
responsible for the dense CSM [15]. 

This conference proceedings summarize the work published in [16], where we perform the 
first systematic search for gamma-ray emission from Type Iln SNe in data recorded by the Large 
Area Telescope (LAT) on-board the Fermi Gamma-ray Space Telescope in the energy range of 
100MeV to 300GeV to probe interaction-powered SNe as cosmic-ray accelerators. We look for 
excess gamma-ray emission correlated in space and time with SNe found in optical surveys such 
as the Palomar Transient Factory (PTF) [17, 18]. We use a sample of 147 Type Iln SNe. A full list 
of SNe used in this analysis can be found in [16]. 

To account for theoretical uncertainties on the expected /-ray flux we perform a model inde¬ 
pendent search. 


2. Gamma-ray Data Analysis 

We analyze 57 month of FerniZ-LAT Pass 7 Reprocessed data. We select events of the source 
class and use the P7REP_SOURCE_V15 instrument response functions. Eirst, we search for 
gamma rays from individual SNe, applying a binned likelihood analysis to a one year time window 
starting 4 weeks before optical peak time. We model the SN spectrum as a power-law leaving 
the normalization and index free. Our background model includes Galactic diffuse and isotropic 
gamma-ray emission in addition to point sources from the 2EGE source catalog [19]. Eor each SNe 
we calculate the test statistic (TS) defined as TS = — 2Alog.if, where A log .if is the difference in 
likelihood between a model without including a source at the SN position (background-only model) 
and a model including a source with a power-law spectrum at the SN position. None of the SNe 
positions yields a significant TS value. In Eig. 2 we compare the SNe TS values with TS values 
obtained from repeating the analysis in random positions in the sky, which allows us to estimate the 
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Figure 1: Schematic picture of a SN exploding in dense circumstellar medium (taken from [2]). 



Figure 2: Left; Distribution of SNe and random position TS values. Right; Cumulative distribution of 
random position TS values. The largest SN TS value is shown as vertical solid red line, it corresponds to a 
p-value of 0.0065 and is compared to the 2(7 and 3(7 value. Considering trials factors, this p-value increases 
to 0.6. 


significance of each SN TS value. None exceeds the 3a level. The largest found TS corresponds 
to a p-value of 0.0065, which increases to 0.6 considering trials factors. 

In a second step, we drop the assumption of a power-law spectral shape for the SNe and instead 
apply a bin-by-bin likelihood analysis, where we allow the normalization of the SN flux to float 
in each energy bin while the background parameter (i.e., background point sources, diffuse and 
isotropic emission) are fixed to values obtained from the broad-band fit. We repeat this analysis 
for three different time windows: 3, 6 and 12 months. The results for one source, SN2010jl, are 
displayed in form of the “Castro plot” in Fig. 3 (similar plots for the other SNe can be found in [16]). 
The Castro plot shows the bin-by-bin likelihood function used to test for a putative gamma-ray 
source at the position of the SN. It allows us to fit any kind of spectral model to the results without 
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Figure 3: Castro plot for SN 2010jl: Within each bin, the color scale denotes the variation of the logarithm 
of the likelihood with respect to the best-fit value of the SN flux using a 12 months time window. 95% CL 
upper limits are shown as black arrows for the 12 months time window and as dotted-dashed and dotted lines 
for a 6 and 3 months time window respectively. For the particular case of SN 2010jl we repeated the analysis 
for an extended time window spanning 4.5 years. The results are overlaid as cyan dotted line. 


repeating the gamma-ray data analysis, i.e. the results are independent of any kind of spectral 
model assumption. The bin-by-bin likelihood results can be used to re-create a global likelihood 
for a given signal spectrum by tying the signal parameters across the energy bins (see [16] and [20] 
for details). 

To be more sensitive to a weak gamma-ray signal we combine the sources in a joint likelihood 
analysis. For this analysis we select the 16 optically brightest and closest sources, fulfilling z < 
0.015 or m < 16.5, where z is the redshift and m the optical peak magnitude^ We repeat the 
joint likelihood analysis for 3 different time windows (3, 6 and 12 months) and for two different 
weighting schemes. In one case we weight each source with l/cf', i.e., we assume that all sources 
have a similar gamma-ray luminosity. In the other case we weight with a weight proportional to 
jQ-o.4m^ i.e., assuming that the gamma-ray luminosity is correlated with the optical luminosity. We 
produce the Castro plot for both cases and fit the results with a power-law and power-law with 
exponential cut-off spectral model (see Fig. 4). None significantly improves the likelihood of the 
fit. 

3. Results 

We use our flux limits presented above to set limits on the gamma-ray luminosity and the 
ratio of gamma-ray to optical luminosity (for details of the limit calculation see [16]). We use a 
gamma-ray emission model based on [1], which assumes that gamma-ray emission is produced 
by CRs accelerated at the early collisionless shock between SN ejecta and the dense CSM. The 
gamma-ray emission can be predicted following the procedure outlined in [2], if the model param¬ 
eters are determined by optical and X-ray observations. However, such model-dependent analyses 

*For some SNe the optical peak magnitude was not determined. In those cases we use the discovery time and 
magnitude. 
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Figure 4: Castro plot for joint likelihood results. Left using a I/if' weighting. Right using a 
weighting. Within each bin, the color scale denotes the variation of the logarithm of the likelihood with 
respect to the best-ht value of the SN flux using a 12 months time window. 95% CL upper limits are shown 
as black arrows for the 12 months time window and as dotted-dashed and dotted lines for a 6 and 3 months 
time window respectively. 


was deferred to future work. Instead in [16] we take a model-independent approach aiming to con¬ 
strain the gamma-ray luminosity as a function of the proton spectral index F. We assume that the 
CR proton spectrum follows a power-law function in momentum and in the calorimetric limit the 
gamma-ray spectral index follows the proton spectral index. We do not take into account the effect 
of gamma-ray absorption, but [2] showed that GeV gamma-rays can escape without severe matter 
attenuation if the shock velocity is high enough. 

As an example we show the predicted gamma-ray flux for the optically brightest source of 
the sample, SN2010jl using a generic gamma-ray flux model. A more detailed modeling of this 
source based on multi-wavelength observations will follow in future work. Here we assumed a 
proton spectral index of F = —2 and a normalization of the gamma-ray flux that yields to a ratio 
of gamma-ray to optical luminosity between 0.01 and 0.1. Furthermore we assumed a distance of 
48.7 Mpc and an apparent R-band peak magnitude of 13.2. Our resulting flux upper limit touches 
the optimistic model prediction, i.e. brings us into reach of the interesting parameter space. Better 
constraints on the gamma-ray escape fraction are crucial to calculate stringent limits on the proton 
acceleration efficiency and will be obtained in more detailed modeling. 

In Fig. 6 (left panel) we show our limit on the gamma-ray luminosity assuming that all sources 
have the same gamma-ray luminosity Ly compared to the limit obtained from the analysis of the 
closest SN in our sample (SN2011ht with a distance of d = 17.7 Mpc). For some choices of F 
the single source limit is better than the combined limit, indicating a statistical under-fluctuation 
in the individual analysis of this source or an over-fluctuation in one of the sources included in the 
joint likelihood. We also show the limit on the ratio of gamma-ray to optical luminosity Lj/Lr, 
where Lr is the R-band luminosity (see Fig. 6, right panel). Note that Lj/Lr is an upper bound on 
Ly/L^^^, where is the bolometric radiation luminosity. We compare the joint likelihood limit 
to the limit obtained using only the optically brightest source of our sample, SN 2010jl. 

4. Conclusion 

The origin of the multi-wavelength emission of Type Iln SNe and the onset of CR production 
in SN remnants is not fully understood. Type Iln SNe are expected to be host sites of paiticle 
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Figure 5: Predicted gamma-ray flux using a generic gamma-ray emission model for SN2010jl including 
theoretical uncertainties (shown as green band) compared to 95% CL flux upper limits obtained from ana¬ 
lyzing a one year time window in blue. 



rp 
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Figure 6: Left: Limits on the gamma-ray luminosity using the joint likelihood results assuming that all 
sources have the same gamma-ray luminosity compared to the limit obtained from analyzing only the closest 
source in our sample (SN2011ht). Right: Limits on the ratio of gamma-ray to optical luminosity based 
on the joint likelihood results assuming a correlation of optical to gamma-ray luminosity compared to the 
single-source limit from the optically brightest source of the sample (SN 2010jl). 


acceleration, which could be pinpointed by transient gamma-ray signals. In this conference pro¬ 
ceedings we presented the results of [16], where we for the first time apply a systematic search for 
gamma-ray emission from a large ensemble of Type Iln SNe in coincidence with optical signals. 
We did not find evidence for a signal, but our limits start to reach the expected model parameter 
ranges. 

We use a joint likelihood analysis to set stringent limits on the gamma-ray luminosity and the 
ratio of gamma-ray to optical luminosity. For proton spectral indices of F < 2.7 we can exclude at 
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95% CL that Lj/Lr >0.1 assuming that Lj/Lr is constant. Future model-dependent calculations 
based on multi-wavelength observations will allow us to also set stringent limits on the proton 
acceleration efficiency. 

The limits presented here are based on minimal assumptions about the gamma-ray production 
and can be used to test various models. 
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